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We consider meta materials made from a two-dimensional dielectric rod-type photonic 
crystal. The magnetic response is studied within the recently developed homogenization 
theory and we in particular study the effects of topology and inhomogeneous broaden- 
ing. While topology itself mainly affects the Mie resonance frequency we find that the 
dispersion in the topological radius S% of the dielectric rods may lead to significant inho- 
mogeneous broadening and suppression of the negative-/! phenomena for 5& ^> e" /e', 
with e = e' + is" being the relative dielectric function of the rods. 
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1 Introduction 

The development of meta materials has led to renewed interest in the original work by 
Veselago on negative refraction [1] and the field of meta materials is in rapid growth. For 
more details we refer to recent reviews [2-4] as well as to recent special issues [5,6]. In 
this work we consider a photonic crystal composed of high-index dielectric rods arranged 
in a two-dimensional lattice, see Fig. ^A). For temporal harmonic modes the H field is 
governed by [7, 8] 

V x _Lv x H(r) = »(r)^H(r) (1) 

e(r) c z 

where e is the relative dielectric permittivity, /i is the relative magnetic permeability, c is 
the speed of light in vacuum, and u = ck is the angular frequency with k = 2tt/\ being 
the free-space wave vector. For Nature's own dielectric materials we locally have, almost 
by definition, that = 1 so that everywhere in space B(r) = fi(r)H(r) = H(r) [8]. 
However, for artificial composite dielectric structures in the meta material limit, where the 
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free-space wavelength A is much longer than the characteristic length scale (periodicity) d 
of the underlying (quasi) periodic structure, one may for the slowly varying homogenized 
fields (envelopes) have an effective permeability fi e g which is constant in space and which 
at certain frequencies deviates remarkably from the microscopic permeability, i.e. /x e fr ^ 
and in particular it may become negative. The same apparent paradox is also possible 
for the dielectric response. As an example Fig. El shows a plane-wave simulation [9] of 
Eq. {1} for the structure in Fig. ^A. with e = 200 and fi = 1 for the rods and e = 
H = 1 for the background air region. For these inverse photonic crystal structures, the 
band structure and the photonic band-gap formation can be understood in terms of Mie 
resonances [10, 11]. Clearly, the band structure reveals a strong hybridization at the 
frequency corresponding to the first Mie resonance of the rods as indicated by the dashed 
line. It is this resonance behaviour which make the rods behave as artificial magnetic 
atoms with the possibility for a negative-/! response. 

In this work we extract the effective permeability by applying the recently derived meso- 
scopic homogenization theory for p-polarized fields by Felbacq and Bouchitte [12-14]. We 
in detail study the effective permeability /i e g and the effect of the topology of the rods as 
well as the influence of spatial disorder/randomness and dispersion in the rod topology, 
compare panels (B) and (C) in Fig. ^ For rods with a more arbitrary cross section Q we 
define a topological radius by 



where &/ is the area of Q and & is the perimeter of its surface dQ, see Fig. [D^B). It 
is easily verified that a circular cross section of radius R gives ffl = R. More generally, 
the topological radius S% captures the qualitative resonance behaviour of the effective 
permeability as we shall see in the following. 

Our results for more arbitrary topologies rely on approximative results for the spectral 
problem of the Helmholtz equation as employed recently by one of us in the context of 
diffusion dynamics, fluid flow, and ion transport in microfluidic problems [15-17]. 

2 Multi-scale expansion theory 

We consider the problem of p-polarized fields corresponding to a magnetic field aligned 
parallel to the rods, see Fig. H^A). In the limit where A ^> d the problem posed by 
Maxwell's equations can be treated by a multi-scale expansion approach leading to an 
effective magnetic permeability given by [13] 



2s// & 



(2) 



/u c ff = / x /i^ + (1 - /) x /in 




(3) 
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where = is the area of f2, / is the air-filling fraction, and m) is governed by 

V 2 |m) + k 2 e\m) = 0, re!) (4a) 
\m) = 1, re dfi. (4b) 

Above, we have used the Dirac bra-ket notation and inner products are defined as 

(h\g) = / drhg, (5) 



where h(r) and g(r) are functions on f2. Similarly, m), or m(r), in Eqs. (3) and (4) 
corresponds to the magnetic field component along the rods and |l) denotes the unity 
function. 

In the following we assume a quadratic unit cell with area d 2 , see Fig. E^B), so that 
/ = 1 — /d 2 . As in Refs. [12-14] we solve the problem in a complete basis defined by 
the Helmholtz eigenvalue equation 

V 2 |^) + ^n) = 0, reQ, (6a) 
\^ n ) = 0, re dn, (6b) 

with the basis functions being orthonormal, i.e. (jj) m \ip n ) = b nm - Obviously, the Ansatz 

H = I 1 ) + ^«n|^n) (7) 

fulfills the boundary condition and the expansion coefficients are readily determined by 
standard bra-ket manipulations. In this way we arrive at the following general expression 
for the effective permeability [12] 

„ -i i v - fc2g kwi 2 r <A 

n n 

The strength of this result is obviously the connection of /x e fr to the spectrum and eigen- 
states of the Helmholtz equation which have been studied in great detail in various con- 
texts of physics, including membrane dynamics, the acoustics of drums, the single-particle 
eigenstates of 2D quantum dots, not to forget quantized conductance of quantum wires. 
Finally, the results of course in turn provides the direct link between /x e fr and the Mie 
resonances. 

3 Approximate properties of the spectral problem 

In the following we utilize approximate results for the spectral problem and its dependence 
on geometry. We start by introducing the effective area s^ n = |(l|-?/> n )| 2 so that Eq. (jHJ) 
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is rewritten as 

-k 2 e g/ n 



/^cfT 



n 



+ (9) 

where we in the second line have truncated the rapidly converging sum making the ap- 
proximation valid in the vicinity of the main resonance. 

In Table we have listed the first eigenvalue k\ and the corresponding effective area 
srf\ for a number of different geometries. For the numerical results we have employed 
a finite-element method with an adaptive mesh algorithm (Comsol multiphysics) ensur- 
ing fast convergence of both eigenvalues and eigenstates. Inspecting the column for the 
eigenvalue k\ in Tabled we note that the numbers do not in general change much from 
geometry to geometry, despite the variation in topology. The small spread originates in 
the normalization by the topological radius ffl, see Eq. (J2J). For the effective area sd\ we 
likewise observe a modest variation. It is worthwhile emphasizing the results for a circular 
cross section which become good approximations for also other relatively compact topolo- 
gies. In the following we thus use the approximation that k\ ~ Xoo/& and srf\j srf ~ 4/xq 
where Xoo — 2.405 is the first zero of the Bessel function, i.e. Jo(Xoo) = 0. With these 
approximations we finally arrive at 

4 -{mfe 



Meff ~ 1 + (1 - /) 



X 2 00 {k^fe - x 2 00 



= l + (l-/)4 -o 2 (10) 

*oos-(^) 2 (A) 2 

where all the dependence on geometry is captured by the topological radius 

We may now estimate the resonance wavelength. Suppose that e = e' + ie", then for 
weakly absorbing materials, e" <C s', we have the resonance around 

27T / — 2tt i — _ . 

A ~ — yfe' « — Ve 1 & (11) 
«i Xoo 

with a line width of the order 

8\~—±-m. (12) 



Xoo Ve' 

As expected the line width reflects the damping due to e" in the high-index material. 
Furthermore, we may also estimate the magnitude of the effective magnetic response. 
Just below the resonance the real part is bounded by 



while for the imaginary part 



ReW^l4(l-/)J (13) 

A00 fc 



HM<4-(i-/)^. (14) 

A00 fc 
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The above results illustrate how e, /, and St influence the magnetic response. Obviously, 
a negative permeability requires that e' ^> e". On the other hand, some absorption is of 
course needed to have a finite line width instead of a delta-function-like response which 
would suffer tremendously from the inevitable effects of inhomogeneous broadening. We 
shall return to this in Sec. El 

4 Comparison to numerical results 

The validity of the assumptions underlying Eq. (JHJ) has been tested in Refs. [12-14] by 
comparing simulated transmission spectra of the real photonic crystal to spectra based on 
the corresponding homogenized meta material. In this work we thus assume Eq. (J3J) valid 
and focus on checking the quality of the approximate result in Eq. (|10|). We do this by 
comparing the outcome of Eq. (jlOJ) to numerical solutions of the original problem posed 
by Eq. (jH)). We use a finite-element method with an adaptive mesh algorithm (Comsol 
multiphysics) for this purpose. For the dielectric function we use e = 200 + 5z as in 
Ref. [12-14]. Fig. El illustrates the corresponding photonic band structure. Clearly, there 
is a strong hybridization at the frequency corresponding to the isolated Mie resonance of 
the rod, see Eq. (Jll|) . It is this avoided crossing which causes strong resonance in the 
effective permeability that is observed in Fig. El The approximate result in Eq. fTT)|) works 
excellently near the main resonance for rods with circular cross section as seen in Fig. El 
while the approximation of course doesn't capture high-order modes due to the truncation 
of the summation in Eq. (JHJ). In the following we go a step further and test the result 
for highly non-circular geometries. Fig. 0] illustrates results for one such example. As 
seen there is a good qualitative agreement and since the geometry is highly non-circular, 
see inset, the quantitative agreement is perhaps better than expected. Note how the 
change in topological radius & scales the resonance wavelength A* and the line width S A 
compared to the circular case, Fig. El in accordance with Eq. (II lj) while the conservation 
of the air-filling fraction / leads to the same magnitude of the response in accordance 
with Eqs. (fT3|) and (|T4^) which depend only on /, but not on 

5 Inhomogeneous broadening 

The present results suggest that shape plays only a minor role and mainly affects the 
particular value of the resonance frequency (which is proportional to the topological radius 
M) while the air-filling fraction / and the dielectric function e = e' + ie" tend to govern the 
magnitude of the effective magnetic response. This suggests that homogenized materials 
will not be strongly influenced by spatial disorder /randomness conserving the average air- 
filling fraction while shape dispersion (dispersion in the topological radius, possibly also 
influencing the average air-filling fraction) will give rise to a smearing of the resonance 
phenomena. Suppose that for a particular topological radius the resonance wavelength 



5 



is Aq = a&o, see Eq. (ITT]) , and that the width of the resonance is SXq. Assuming a Gaussian 
distribution of 



oc exp 



(15) 



and neglecting the dependence of the intrinsic line width 5Xq on 3& we get an inhomoge- 
neous suppression of the resonance magnitude by a factor II, 



n 



(16) 



and an inhomogeneous broadening of the resonance by a factor II 

5X ~ II~ 1 (5Ao. 



(17) 



Thus, inhomogeneous broadening will be negligible if a ( ) <C 1 corresponding to 



W « T (18) 

Fig. El illustrates these results in the case of rods of circular cross section as in Fig. El 
The results are obtained from a numerical evaluation of J °° ' dM \i e ^{J%)P{S^) with /x e g 
given by Eq. (jlOj) . As also suggest by Eq. (|TH|) a size dispersion with a relative variation 
comparable to e" / e' preserves the negative-/! phenomena while larger variations will result 
in a positive /x e g. 



6 Conclusion 

We have studied the magnetic response of rod-type dielectric meta materials within the 
recently developed homogenization theory [12] and we have in particular investigated the 
effects of topology and inhomogeneous broadening. We predict that topology itself mainly 
affects the Mie resonance frequency while the resonance shape mainly depends on the air- 
filling fraction and the ratio e'/e" with e = e' + is" being the relative dielectric function 
of the rods. Furthermore, we argue that randomness/disorder conserving the air-filling 
fraction will have a limited influence on the negative-/^ phenomena while dispersion in the 
topological radius & of the rods may lead to significant inhomogeneous broadening and 
suppression of the phenomena. 
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Table 1: Central dimensionless parameters for different geometries. a See e.g. [18] for 
the eigenmodes and eigenspectrum. b Data obtained by finite-element simulations. c See 
e.g. [19] for the eigenmodes and eigenspectrum. 
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Figure 1: (A) Photonic crystal consisting of parallel high-index dielectric rods. We con- 
sider the p-polarized case with the magnetic field parallel to the rods. (B) Top view of 
the photonic crystal with rods arranged in a periodic lattice with pitch d. (C) Top view 
of weakly disordered/random photonic crystal. 
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r x m r 

Figure 2: Plane- wave simulations of photonic band structure for a d x d square lattice 
of rods of circular cross section with e = 200, see Fig. E^B). The air-filling fraction 
/ ~ 80.37% corresponding to a rod of radius R — & — d/4. The dashed line shows the 
fundamental Mie resonance of an isolated dielectric rod, Eq. fTTf . 
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Figure 3: Effective permeability /i c fj versus normalized wavelength X/d for a rod with 
e = 200 + 5i in a quadratic d x d unit cell with an air-filling fraction / ~ 80.37% 
corresponding to a rod of radius R = M = d/4. The inset illustrates the unit cell as well 
as the function m at a wavelength above resonance. The solid lines are numerical exact 
solutions based on a finite-element (adaptive mesh) solution of Eq. (j3J) while the dashed 
lines show the approximate result in Eq. (jlOj) . 
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Figure 4: Effective permeability /i c fj versus normalized wavelength X/d for a rod with 
e = 200 + 5i in a quadratic d x d unit cell with an air-filling fraction / ~ 80.37%. The 
inset illustrates the unit cell as well as the function m at a wavelength above resonance. 
The highly non-circular cross section has a topological radius of 3& ~ 0.151 x d. The solid 
lines are numerical exact solutions based on a finite-element (adaptive mesh) solution of 
Eq. (J3J) while the dashed lines show the approximate result in Eq. The horizontal 

dashed lines show the limits in Eqs. (TT3J) and (JHJ). 
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Figure 5: Inhomogeneous broadening of the effective permeability /x e g for rods with a 
circular cross section in a quadratic d x d unit cell with an average air-filling fraction 
/o ~ 80.37% corresponding to rods with an average radius R = & = d/A. The dashed 
curves show the intrinsic line width corresponding to the results in Fig. El 
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